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ABSTRACT 



Context. Subcritical transition to turbulence in Keplerian accretion disks is still a controversial issue and some theoretical progress 
is required in order to determine whether or not this scenario provides a plausible explanation for the origin of angular momentum 
transport in non-magnetized accretion disks. 

Aims. Motivated by the recent discoveries of exact nonlinear steady self-sustaining solutions in linearly stable non-rotating shear 
flows, we attempt to compute similar solutions in Rayleigh-stable rotating plane Couette flows and to identify transition mechanisms 
in such flows. 

Methods. Numerical simulations (time-stepping approaches) have been employed in most recent studies of the problem. Here, we 
instead combine nonlinear continuation methods and asymptotic theory to study the nonlinear dynamics of plane Couette flows in 
Rayleigh-stable regimes. 

Results. We obtain exact nonlinear solutions for Rayleigh-stable cyclonic regimes using continuation techniques proposed by Waleffe 
and Nagata but show that it is not possible to compute solutions for Rayleigh-stable anticyclonic regimes, including Keplerian flow, 
using similar techniques. We also present asymptotic descriptions of these various problems at large Reynolds numbers that provide 
some insight into the differences between the non-rotating and Rayleigh-stable anticyclonic regimes and derive some necessary con- 
ditions for mechanisms analogous to the non-rotating self-sustaining process to be present in (constant angular momentum) flows on 
the Rayleigh line. 

Conclusions. Our results demonstrate that subcritical transition mechanisms cannot be identified in wall-bounded Rayleigh-stable an- 
ticyclonic shear flows by transposing directly the phenomenology of subcritical transition in cyclonic and non-rotating wall-bounded 
shear flows. Asymptotic developments, however, leave open the possibility that nonlinear self-sustaining solutions may exist in un- 
bounded or periodic flows on the Rayleigh line. These could serve as a starting point to discover solutions in Rayleigh-stable flows, 
but the nonlinear stability of Keplerian accretion disks remains to be determined. 
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1. Introduction 

Physical mechanisms responsible for turbulent angular momen- 
tum transport in Keplerian accretion disks have long been a 
matter of debate among astrophysicists. While there is con- 
vincing evidence t hat the presence of magnetic fields in disks 
dDonati et al.l2005l) triggers turb ulence thanks to a linear magne- 
torotational instability (MRJ, s ee IVelikhovll959t IChandrasekhaJ 
U960tlBalbus & Hawlevlll991l) . purely hydrodynamic sources of 
turbulence still have to be discovered in some regions of pro- 
toplanetary disks which are like ly to be decoupled from mag- 
netic fields (iFromang et a l. 2002) and in disks forming in binary 
systems in quiescence, in which MHD turbulence and the asso- 
ciated turbulent transport may die away as a consequence of the 
small ness of the magnetic Reynolds number dGammie & Menoul 
fl998h . 

Studying hydrodynamic transition to turbulence in Keplerian 
flows is an extremely difficult task for essentially two rea- 
sons. On the one hand, there is currently no known local 



linear instability which could provide a route to turbulence 
in the non-magnetized Keplerian flow regime. iDubrulle et al.l 
d2005l) have pointed out the existence of an instability 
in stably str atified accretion disks, but its local nature i s 
very unclear dUmurhanl l2006t iBrandenburg & Dintransll2006l) . 
On the other hand, the linear operators associated with 
shear flows are mathematically non-normal, which physi- 
cally means that short-term dynamics (transient algebraic 
growth) rather than normal modes may play a n important 
role in the process of transi ti on to turbulence (Korvcanskyl 
1992h iTrefethen et alJ Il993l: ISchmid & Henningsonl I 2000 ; 
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Chagelishvili et alJ l2003t iTevzadze et alJ l2003t lYeckol I2004 
Afshordi et all l2005t iMukhopadhvav et alJ |2005|) . Several au- 
thors dRi chard & Zahn|[l999t lLon garettH |2002|) have suggested 
that subcritical transition may be possible in accretion disks in 
spite of the stabilizing effect of the Coriolis force , but this point 
of view ha s been conteste d by Balbus et al. ( 1996i). lHawlev et aTl 
d!999t) and lBalbusI d2003l) . The relevance of such mechanisms to 
the proble m of turbulent transport in disks has also been ques- 
tioned by iLesur & LongarettT d2005l) . who extrapolated from 
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numerical simulation results that the critical Reynolds number 
for this kind of process to occur in the Keplerian regime was 
far too large and the associated angular momentum transport 
far too small to accoun t for the accretion process. Recently, 
iBalbus & Hawlevl d2006l) have further questioned this transition 
scenario, basing their argument on the observation that tran- 
siently amplified shearing waves are solutions for any ampli- 
tude (because their self-interactions are vanishing) and canno t 
ther efore direct l y trigg er turbulence. IBalbus & Hawlevl d2006l) 
and IShen et al.l (120061) have also studied the nonlinear evo- 
lution of these non-axisymmetric shearing waves using high- 
resolution simulations and have found that they become unsta- 
ble to Kelvin-Helmholtz instabilities if they reach a sufficient 
amplitude. However, the turbulent motions generated this way 
ultimately die away in the configuration of their simulation. 

Most recent theoretical studies of hydrodynamic turbulence 
in accretion disks have focused on linear transient growth, which 
is only one part of a full subcritical transition scenario, and is not 
a sufficient condition for nonlinear instability. Nonlinear dynam- 
ics has mostly been stu died numerically within the shearing box 
model dRogallol[T98l using direct time-stepping approaches. 
The nonlinear outcome of the transient growth process remains 
unclear and theoretical progress is required in order to clarify the 
existence and nature of subcritical transition in Rayleigh-stable 
(centrifugally stable) rotating shear flows. The aim of this pa- 
per is to present some theoretical and numerical arguments that 
may help to make progress in this direction, and notably to in- 
vestigate if some nonlinear mechanisms that are crucial for sub- 
critical transition to turbulence in non-rotating shear flows play 
a role in linearly stable rotating shear flows too. We have been 
motivated by the vast amount of experimental, numerical and 
theoretical work that has become available in fluid dynamics in 
the last ten years and have oriented our study towards one of 
the most promising avenues of research i n that field, namely the 
self-sustaining process, first identified bv lHamilton et al.l d 1995b 
for non-rotating plane shear flows. Our hope has been to be 
able to identify a similar process for Rayleigh-stable rotating 
shear flows, notably the anticyclonic plane Couette flow on the 
Rayleigh line, which at first glance presented many similari - 
ties with its non-rotating counterpart dBalbus & Hawlevlll998l) . 
Unlike most recent studies, we have tried to address the fully 
nonlinear problem using continuation methods, in order to com- 
pute accurate nonlinear solutions in different rotation regimes 
and to pinpoint the physical mechanisms that may give birth to 
turbulence in these regimes. Our results shed some new light on 
the nonlinear outcome of linear transient growth on the Rayleigh 
line and show that it is not possible to identify nonlinear solu- 
tions in wall-bounded Rayleigh-stable anticyclonic shear flows 
using arguments similar to those pertaining to the self-sustaining 
process in non-rotating shear flows. We also argue that these 
findings do not rule out the possibility of a subcritical transition 
scenario in accretion disks, but show that self-sustaining mech- 
anisms in rotating plane shear flows with linearly stable anticy- 
clonic rotation, if they exist, must be qualitatively different from 
their non-rotating equivalent. 

The outline of the paper is as follows. In Sect. [2] we sum- 
marize the current understanding of transition to turbulence in 
shear flows and give a detailed description of the phenomenol- 
ogy of the self-sustaining process for non-rotating shear flows. 
Section[3]is devoted to the mathematical description of the non- 
linear problem that we have been attempting to solve. The nu- 
merical tools developed for the purpose of our study are pre- 
sented in Sect. [4] We then report some numerical results for sub- 
critical transition in cyclonic (Sect. [5]) and anticyclonic (Sect.[6]l 



rotating plane Couette flows. Section|7]provides some theoretical 
insight into the asymptotics of self-sustaining processes at very 
large Reynolds numbers in order to understand these numerical 
results and the difficulties in identifying a self-sustaining pro- 
cess so far for linearly stable anticyclonic regimes. A discussion 
of the results and guidelines for future work are finally presented 
in Sect. [8] 



2. Transition to turbulence in shear flows and the 
self-sustaining process 

To gain some theoretical understanding of the Keplerian prob- 
lem, it is necessary to come back to the more general problem 
of hydrodynamic transition to turbulence in shear flows. For a 
long time, the only classical non-rotating shear flow for which an 
explanation for transition was available is plane Poiseuille flow 
with rigid boundaries. The transition was attributed to the linear 
inst ability of this flow at a critical Reynolds number Re = 557 2 
(see lOrszaddl97ll) : IZahn e"taildl974l) : lbrszag & Paterald 19801) ). 
which takes the form of travelling waves. The instability is sub- 
critical and finite amplitud e solutions can be tracked down to 
Re = 2900 (lHerbertlll976l) . However, it was noticed that this 
global critical Reynolds number was rather high compared to the 
transition value observed experimentally. Studies of secondary 
instabilities of these t ravelling waves di d not seem to resolve this 
discrepancy (see e. g. lDrissi et alj!999b . Another major problem 
with this approach was that it applied only to plane Poiseuille 
flow: plane Couette flow and pipe Poiseuille flow do not have 
any linear instability and nevertheless become turbulent at mod- 
erate values of Re (the linear stability for all Reynolds numbers 
of Plane Couette flo w with rigid boundaries was demonstrated 
by |Romanovldl973l) ). 

In the last decade, a nonlinear self-sustaining process act- 
ing as an organizing center of the dynamics in phase space has 
been discovered, leading to a new paradigm in the understand- 
ing of transition in shear flows. Corresponding exact nonlin- 
ear solutions to the Navier-Stokes equations have been iden- 
tified numerically in plane Poiseuille flow and plane Couette 
flow with either no-slip or str e ss- free boundaries ( Hamilton et al.l 
1991 IWaleffel fl99l fl997t ISchmiegell 119991: IWaleffel l200ll 



2003b. and in pipe Poi seuille flow (IFaisst & Eckhard 



2003; 



Wedin & Kerswellll2004l) . leading to an elegant solution to the 
historical problem of transition in pipes dRevnoldsll 18831) . These 
coherent structures take the form of nonlinear travelling waves 
or of a steady pattern, depending on the symmetries of the 
original problem. The critical global transition Reynolds num- 
bers and the mean flow profiles for these solutions prove to 
be far closer to experimental measurements than predictions 
made by other theories, and an experimental detection of these 
structures in pip e Poiseuille flow has been reported recently 
dHof et al.ll2004l) . The self-sustaining process as it is understood 
now is a subtle balance between linear and nonlinear mecha- 
nisms. The linear part of the process results physically from 
the transient growth of a streamwise-independent (axisymmet- 
ric in the accretion disk terminology) streamwise velocity field 
(referred to as the streaks field) which originates from the redis- 
tribution of the mean shear by low-amplitude streamwise vor- 
tices (or equivalently streamwise rolls). Such a growth is lim- 
ited to time scales 0(Re) because of the ultimate diffusion of 
streamwise vortices. This mechanism is called the lift-up effect 
dLandahllll980l:lButler & Farrelllll992l) . A steady state can only 
be achieved thanks to the development of a three-dimensional 
linear instability of the high-amplitude streaks, which produces 
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a nonlinear feedback that is sufficient to maintain the orig- 
inal rolls field. Note that the linear growth of streamwise- 
independent structures is a key requirement for the whole pro- 
cess to occur, and that this physical mechanism is a direct con- 
sequence of the non-normality of the linear operator. It is also 
different from and more powerful than the "swing amplifica- 
tion" of stream wise-dependent (non-axisymmetric) Orr-Kelvin 
shearing waves (|Lord Kelvinlll887t IOrrlll907|: | Knoblochlll984t 
ICraik & Criminalell986l:lButler & Farrellll 19921)" It is worth em- 
phasizing here that the upper and lower branches of these three- 
dimensional solutions are both unstable to steady and oscilla- 
tory modes wi th the same basic symmetries. As discussed by 
IWaleffd d2001l) . this does not mean that the original coherent 
structures are not relevant to transition. On the contrary, they 
are expected to show up intermittently in the flow, and the fact 
that they have instabilities ensures that bifurcations to even more 
complex structures and to turbulence are possible. 

Far less is known analytically regarding the stability of span- 
wise rotating plane shear flows (the spanwise direction is the one 
perpendicular to the disk plane in an accretion disk). The linear 
stability properties of such systems depend on the rotation num- 
ber Rq, which is proportional to the ratio of the rotation rate O 
and the background shear rate S (defined here as the negative of 
the background flow relative vorticity): 



Rn = - 



2Q. 



(1) 



This R q is the same as the one defined by iLesur & Longarettil 
d2005l) (S here is —S in their paper but the y use Rn — 2D./S) 
and d iffers in sign from the definition of iBech & Anderssonl 
(11997b . Here we assume that S is constant in the shearwise 
(radial) direction, which corresponds to a plane Couette flow. 
A necessary condition (known as the Rayleigh criterion) for 
this flow to become linearly unstable is that -1 < Rp_ < 0, 
which shows that the flow has to be anticyclonic to be un- 
stable, i. e. that the vorticity of the background flow and the 
rotation rate have opposite sign. In this so-called Rayleigh- 
unstable regime, the centrifugal instability is responsible for 
the transition to turbulence, as observed in both labora tory 
dTillmark & Alfredssonlll996HAlfredsson & Tillmarkll2005l) and 
numerical dBech & Anderssonl 1 19961 1 19971) exp eriments. This 
instability is the same as the instability studied by iTavlorid 1923b 
in the Taylor-Couette experiment, which exhibits axisymmetric 
structures called Taylor vortices close to the bifurcation point. 
On the contrary, rotating plane Couette flow in the cyclonic 
regime (Rn > 0) and anticyclonic flows with Rq < -1 are 
Rayleigh-stable for all Re. Local linearizations of Keplerian 
flows, which have Rn = -4/3, are included into this category. 
From the experimental point of view, turbulence is of course 
known to be present in plane Couette flow for Rn = (non- 
rotating plane Couette flow), but it has also been observed fo r 
small positive rotation numbers dTillmark & Alfredssonl 1 1996b . 
Note finally that unlike for non-rotating plane Couette flow, there 
is no complete proof that rotating plane Couette flow at Rn = - 1 
is linearly stable for all Reynolds numbers. 

At that point, it should be noted that a strong theoretical con- 
nection exists between the nonlinear stability of non-rotating 
plane Couette flow and the weakly nonlinear behaviour of 
Rayleigh-unstable Taylor-C ouette flow in the na rrow gap limit. 
This link was discovered bv lNagataldl988lll990l) . He found that 
the three-dimensional nonlinear self-sustaining steady solutions 
for non-rotating plane Couette flow (which were unknown at 
that time and had therefore not been related to any kind of self- 
sustaining process) could be obtained directly by nonlinear con- 



tinuation (with respect to the rotation rate) of wavy Taylor vor- 
tices, which emerge from a three-dimensional bifurcation of the 
Tayl or vortex flow. For a given range of streamwise wavenum- 
bers. lNagatal ([1990) showed that these solutions could be contin- 
ued slightly into the cyclonic regime. iKomminaho et alj d 19961) 
performed numerical simulations and discovered that the large- 
scale structures of turbulence disappeared in the cyclonic regime 
at Rn = 0.06 for Re - 750, in qualitative agreemen t with the ex- 
perimental results of iTillmark & Alfredssonl dl996l) . 



Whether or not turbulence exists in the Rayleigh-stable, an- 
ticyclonic regime (and notably in Kep lerian accretion disks) i s 
still very unclear, as mentione d earlier . iRichard & Zahnld 1999b . 
using data from experiments bv lWenda d 1 933b and lTavlorld 1936b . 
suggested that transition to turbulence should be observed in 
a Taylor-Couette flow in that regime at rather modest <9(10 5 ) 
Reynolds numb e rs, and claimed exp e riment al evidence for this 
dRichardl 1200 lb . IBech & Anderssonl d 1997b found that turbu- 
lence disappeared at Rn = -1 for Reynolds numbers up to 
Re ~ 1300. Nonlinear s tability in the Taylor-Cou ette problem 
was also investigated by iGaraud & Ogilvie! d2005b using a clo - 
sure model for the Reynolds stresses derived bv lOgilvidl2003b . 
They showed that within the framework of this model, nonlinear 
solutions in the anticyclonic Rayleigh-stable regime should exist 
if solutions in the cyclonic regime exist, thanks to a symmetry 
with respect to Rn = -0.5. Streamwise-independent solutions 
of t he Navier-Stokes e quati ons do have such a symm etry (see 
also lBalbus et al.ll 1996b . but lLesur & Longarettil d2005h pointed 
out that it is broken when fully three-dimensional solutions of 
the Taylor-Couette problem are considered, leaving the issue of 
the existence of nonlinear solu t ions i n the stable anticyclonic 
regime unsettled. lHawlev et a l. (1999) an d ILesur & Longarettil 
d2005b . using shearing box simulations dRogallol 11981 ). ob- 



serv ed turbulence for an ti cyclon ic rotation down to Rn = - 1 03 
and ILesur & Longarettil d2005b . using either shearing box or 
rigid-boundary simulations, found turbulence at small positive 
rotation numbers (up to 0.12). The nature of t he numerical 
code ( which does not include explicit diffusion) of lHawlev et al.l 



eft 

d 1999b did not allow them to make statements on critical 
Reynolds nu mbers for the onset of turbul ence in these flows, but 
the results of ILesur & Longarettil d2005|) show g o od ag reement 
with the experiments of ITillmark & Alfredssonl d 19961) on the 
cyclonic side and seem to indicate that transition occurs on the 
anticyclonic side for Re ~ 1500 at Rn = -1- Also, their results 
show that the critical Reynolds number for the appearance of tur- 
bulence at more negative rotation numbers increases very rapidly 
(Re c ~ 4 x 10 4 for Rn = —1.03), which is not consisten t with 
the experimental results of lRichardl d2001b . |jiet al.1 d2006b have 
recently performed a new Taylor-Couette experiment. Their re- 
sults show that quasi-Keplerian flows remain essentially steady 
at Reynold s numbers up to millions, i n agreement with the con- 
clusions of lLesur & Longarettil d2005h . 



This summary of our current understanding of transition in 
non-rotating and rotating shear flows makes it clear that we still 
do not apprehend the whole process in a satisfying way. Some 
important questions, such as the existence of a some qualitative 
symmetry with respect to Rn = -0.5 or the existence of a self- 
sustaining process in the presence of rotation, remain largely 
unanswered. From our point of view, these issues have to be ad- 
dressed theoretically before conclusions on hydrodynamic tran- 
sition to turbulence in Rayleigh-stable rotating shear flows such 
as Keplerian flows can be made. 
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3. Setting up the problem 

3.1. Physical model 

In order to address this problem, we considered an incompress- 
ible spanwise-rotating plane Couette flow between rigid bound- 
aries, which is representative of shear flows that are linearly sta- 
ble for all Reynolds numbers in the absence of rotation. As is 
common in the fluid dynamics literature, (x, y, z) stand respec- 
tively for the streamwise (azimuthal, -y in the shearing sheet 
coordinates), shearwise (radial, x in the shearing sheet coordi- 
nates) and spanwise (vertical) directions. The shearwise direc- 
tion is bounded by two rigid walls, so that —d < y < d, and 
d is chosen as unit of length. We further assume that the x and 
z direction are periodic with spatial periods L x — 2nd I a and 
L z = 2nd//3 . The background flow reads Ub = Sye x , the rota- 
tion vector is given by £1 — £1 e z and the constant shear rate S 
is used to define a time unit. In nondimensional form we there- 
fore have -1 < y < 1 and Ub — ye x . The classical definition of 
the Reynolds number in plane Couette flow is based on the shear 
velocity at the boundaries and half of the channel width: 



Re = 



\S\d 2 



(2) 



where v is the constant kinematic viscosity of the fluid. The 
Reynolds nu mber defined a bov e is four times smal l er tha n the 
ones used bv lNagatal d 1986b and lLesur & Longarettil d2005l) . and 
the rotation number, defined in Eq. (Q]), obeys Rq = -f2xc/R e TC 
where the TC inde x refers to qua ntities defined in the Taylor- 
Couette problem bv lNagatal d 19861) . 

3.2. Equations 

The full velocity field, denoted by U, is split into U b, the orig- 
inal streamwise background flow, and velocity perturbations 
u = (u, v, w), which include any possible modification of the 
original background flow. The nondimensional equations gov- 
erning the evolution of U and of the total pressure P = Pb + p 
(implicitly divided by the constant density p), where p denotes 
the pressure perturbation, are the incompressible Navier-Stokes 
equations, which read 



— + U VU -R n e-xU = -VP + —AC/, 
dt z Re 



V 17 



0. 



(3) 
(4) 



These equations have to be completed by periodic boundary con- 
ditions in x and z and by rigid (no-slip) boundary conditions in 
the y direction (the x, z dependence is omitted here): 

u(-l) = u(l) = v(-l) = v(l) = w(-l) = w(l) = 0, (5) 

Anticipating the numerics, an important point is that these equa- 
tions reduce to two three-dimensional scalar equations govern- 
ing the evo lution of shearwise velocity ( v) and vorticity (77) per- 
turbations ( ISchmid & Henningsonl2000h : 

(jt ~ ri a ) Av + Pv {U vv ~ Rnez xU) =°> (6) 

[jt ~ r\ A ) V + p v(V-VU- Rn e z X V) = 0, (7) 



where 



P Y = -e v -VxVx(0, 
P„ = e y ■ V x (•), 



(8) 
(9) 



plus two streamwise and spanwise-independent equations de- 
scribing the evolution of the mean streamwise and spanwise ve- 
locity fields: 



d 


1 d 2 




Redy 2 


d 


1 d 2 


Jt 


Redy 2 



(w) + (e z ■ U ■ VZ7) = 0, 



(10) 



(11) 



where brackets denote (x,z) averages. Note that (v) = re- 
sults from the incompressibility constraint and boundary con- 
ditions ©, and that the Coriolis force does not enter Eqs. ( fTOl - 

Cld. 



3.3. Symmetries 

The fluid motions that will be studied here have several symme- 
tries, which can be used to reduce the computational costs. We 
will first deal with three-component two-dimensional (stream- 
wise independent) flows (hereinafter 2D-3C), like combined 
streamwise rolls (vortices) and streaks, which have a reflection 
symmetry with respect to the z = plane for a properly chosen 
phase in z. This symmetry reads 



Z : (x,y,z) -> (x,y,-z), (y,rf) -> (v,-tj). 



(12) 



Thanks to the X symmetry of the basic 2D-3C flows, fully three- 
dimensional nonlinear solutions bifurcating from these flows 
have either the X symmetry or the so-called "shift-and-reflect" 
symmetry given by 

S 2 :{x,y,z)^{x + L x l2,y,-z), {v,rj) -» (v.-ifl. (13) 

A basic 2D-3C solution has both Z. and S2 symmetries but the 
three-dimensional part of the total velocity field has only one 
of them. Finally, choosing the x-phase of the three-dimensional 
solution (brackets denote (x, y) averages) such that 



3 (^exp(-i'a'x)) = or 3 (vexp(-/cwc)) = 0, 



(14) 



depending on whether the solution has S2 or Z, symmetry, it is 
possible to show that for plane Couette flow the velocity field 
has the following "shift-and-rotate" symmetry: 

Si : (x,y,z) -» (LJ2-x,-y,z+L z /2), (v,i/) -» (-v,i?). (15) 

The forementioned sym metries have been named according to 
the nomenclature used in lWedin & KersweH (|2004|) with (x, y, z) 
here corresponding to (z, r, <p) in their paper on pipe Poiseuille 
flow (not to be confused with the cylindrical coordinates some- 
times used for accretion disks). 



3.4. Comments on the model 

Before we start describing our results, let us briefly comment on 
our choices regarding this physical model. We will discuss two 
points: incompressibility and boundary conditions. Accretion 
disks are fully compressible and our assumption of incompress- 
ibility may appear as a limitation. We fully support the idea 
that sound waves may play an important role in the nonlinear 
dynamics of accre tion disks, notably by destabilizing coheren t 
vortical structures dBroadbent & Moorelll979l:lBodo et alj|2005l) 
and by limiting their extent because of the need to maintain 
causal acoustic contact, and that compressible physics should 
be included in any accretion disk study aiming at being real- 
istic. However, the fact that disks are compressible does not 
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mean that incompressible dynamics is unimportant in this envi- 
ronment. Our motivation in this study has been to identify such 
fundamental incompressible processes, which are major actors 
of the transition to turbulence in non-rotating shear flows. The 
self-sustaining process described in Sect. [2]does not require any 
compressible effect and is already fairly complex, so our choice 
represents a natural simplification to tackle the corresponding 
rotating problem. 

The second point has to do with our choice of a wall-bounded 
flow model, which may not be relevant to the physics of ac- 
cretion disks, notably because walls are responsible for the ap- 
pearance of boundary layers and global modes. There are two 
reasons for this choice here: first, investigating exact nonlinear 
dynamics using continuation methods usually requires that so- 
lutions be steady in a given frame of reference, and finite shear- 
ing boxes do not admit steady non-axisymmetric (x-dependent) 
solutions (non-axisymmetry is mandatory to have sustained hy- 
drodynamic turbulent angular momentum transport in Rayleigh- 
stable and Rayleigh-neutral flows). In fact, using our set of co- 
ordinates and dimensional units, the shearing-periodic boundary 
condition on u(x, y, z, t) is 



u(x,y + Ly,z, t) = u(x - SL y t,y,z, t). 



(16) 



If u(x,y,z, t) depends on x (its first argument) but not on t (its 
fourth argument) then Eq. dT&b cannot be satisfied because the 
r.h.s. depends on time, while the l.h.s. does not. For this reason, 
the finite shearing box does not offer a practical framework to 
investigate the exact nonlinear dynamics of shear flows. It may 
admit non-axisymmetric solutions that are periodic in time, but it 
seems technically difficult to obtain them using standard numer- 
ical methods. The second reason has to do with the phenomenol- 
ogy of rotating shear flow turbulence in laboratory. At first 
glance, the fact that only wall-bounded flows admit steady three- 
dimensional nonlinear solutions may be problematic for self- 
sustaining process theories in accretion disks. However, we note 
that the phenomenology of the onset of turbulence looks very 
similar in unbounded and wall-bounded flows: for instance, tran- 
sition Reynolds numbers observed for cyclonic rotation in shear- 
ing box es or in configurations wi th rigid boundaries are com- 
parable dLesur & LongarettHl2005l) . Also, experiments on wall- 
bounded Rayleigh-stable anticyclonic plane shear flows face ex- 
actly the same difficulty as shearin g box simulations, i. e. turbu- 
lence dies away for Rg cl ose to - 1 ( Tillmark & Alfredssonl 19961 
iBech & Anderssonll 19971) at Reynolds numbers (9(10 3 ) for both 
types of set-ups. Studying the analytically and numerically sim- 
pler case of wall-bounded flows may thus prove useful to under- 
stand the unbounded problem too. 

4. Numerical methods 

Bearing in mind the symmetry of Couette flow with respect 
to the centre of the channel, we aimed at computing steady 
nonlinear solutions of the previous set of equations (in plane 
Poiseuille flow, which does not have this symmetry, solu- 
tions would be travelling). For this purpose, we developed two 
codes: a full n onli near continuation code sim ilar to those of 
IWaleffd d2003l) and IWedin & Kerswelll d2004 to compute ex- 
act three-dimensional nonlinear steady solutions of Eqs. ©- 
(f6T> - <T71> - (fTOb - dTTb thanks to a Newton iteration algorithm, and 
a two-dimensional linear stability code to find accurate ini- 
tial guesses for the nonlinear code. In both codes, discretiza- 
tion is done on a Gauss-Lobatto grid (Chebyshev collocation) 
in the shearwise direction and a Fourier spectral representa- 
tion is used in the x and z directions. Boundary conditions are 



implemented usin g differentiation matrices co mputed with the 
DMSuite package dWeidemann & Redd v 2000). In the continu- 
ation code, nonlinear terms are computed in physical space and 
the steady nonlinear equations are solved in spectral-physical- 
spectral space using a direct LAPACK solver to perform Newton 
iterations. The stability code makes uses of either the direct 
LAPACK eigenvalue solver, which computes the whole spec- 
trum of eigenvalues of the discretized operator, or the iterative 
Arnoldi ARPACK solver to compute individual eigenmodes. 
The nonlinear code was first checked with a two-dimensional 
subcritical test problem, namely the Qrr- S ommerf eld problem 
for a no-slip plane Poiseuille flow dZahn et al.l Il974t iHerbertl 
1 19761) and the full three-dimensional version was successfully 
tested by reproducing the bifurcation diag ram for no-slip non- 
rotating plane Couette flow obtained by Waleffe (2003j) for 
(a,fi) = (0.577, 1.15) (see below). Dealiasing was implemented 
but did not prove necessary to obtain accurate solutions in most 
cases. The linear solver was tested by calculating the limits and 
critical wavenumbers for the centrifugal instability and by quan- 
titatively reprod ucing results obtained for the inflection al insta- 
bility of streaks dWaleffdl 19951 1 1 9971: iReddv et al.ll 19981) in both 
plane Couette flow and plane Poiseuille flow, for which the in- 
stability of streaks takes the form of a travelling wave (non- 
zero eigenfrequency). Pseudo-arclength continuation was imple- 
mented to follow branches of solutions. The previously men- 
tioned symmetries can be turned on or off in the nonlinear code 
in order to reach higher resolutions. We checked that the intro- 
duction of the symmetries module of the code gave the same re- 
sults as the full problem for mild resolutions. We would finally 
like to emphasize here that the numerical methods used in our 
calculations for spatial differentiation have spectral (exponen- 
tial) convergence, and that the maximum affordable resolutions 
(N x , N y , N z ) = (32, 24, 32) used in this paper (which require ap- 
proximately one Gigabyte of memory) allow for the computation 
of fairly complex velocity fields with an excellent precision. 



5. Cyclonic rotation 

5.1. Computation of exact nonlinear solutions in non-rotating 
plane Couette flow 

In this section, we report results obtained for cyclonic rotation. 
In order to reach this regime, we first computed steady non- 
ro tating plane Co uette flo w solutions using the same technique 
as IWaleffd d2003l) (see also lWedin & Kerswelll d2004 for a sim- 
ilar approach for the pipe Poiseuille flow problem). The maxi- 
mum resolution used for these calculations wa s (N x ,N y ,N 7 ) = 
(32, 24, 32), which is similar to the resolution of IWaleffd d2003l) . 
and the Si and £2 symmetries were used in the code to compute 
fully three-dimensional solutions. 

The basic idea of this forcing approach, as we will call it in 
the remainder of the paper, is to compute fully three-dimensional 
nonlinear solutions of the forced Navier-Stokes equations using 
Newton iteration, then decreasing progressively the forcing term 
until a solution for zero forcing is found. In practice, an ini- 
tial guess is required to perform this calculation. In this study, 
we chose (or, /J) = (0.577, 1.15), for which the global critical 
Reynolds number for the self-sustaining process is known to be 
minimized, and set Re = 150. The determination of the ini- 
tial guess was done in two steps. The first step was to compute 
the 2D-3C streamwise-independent part of the flow, denoted by 
(w2D, V2d, W2d) and the second step was to compute the three- 
dimensional instability of this 2D-3C flow. We first focus on the 
first step, which relies on the physics of the lift-up effect in or- 
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der to calculate a combined rolls and streaks velocity field. The 
simplest streamwise rolls structure is given by the most slowly 
decaying solution of the linear Stokes problem for the stream- 
wise vorticity a> = exp (At /Re) a> s (y, z) 



A a> s — A w s 



(17) 



where A is negative. For rigid no-slip boundary conditions and a 
z-periodicity of the rolls [5 = 1.15, the least stable eigenvalue is 
given by A — -9.27306. The corresponding Stokes velocity field 
is 



(Vj,w s ) = (-d z A u> s ,d y A cj s ). 



(18) 



The physics of the lift-up effect is such that (9(1 /Re) streamwise 
rolls generate 0(1) streamwise velocity, so we chose the ampli- 
tude of the rolls to be 



(V2D, W 2D ) = 



A (v s ,w s ) 
Re max(v. s ) ' 



(19) 



where A is assumed to be 0(1). As shown by the sign of A, this 
rolls field decays viscously, so that the following 0(1 /Re 2 ) forc- 
ing term F 

(Fx, F y , F z ) = — (0, v 2D , w 2D ) = —5- t\ ( 2 °) 

Re R e - max(v s ) 

had to be included in the Navier-Stokes equations in order to 
compensate for viscous diffusion. The determination of the 2D- 
3C part of the initial guess was finally completed by solving the 
linear problem 



1 d d\ du B 

— A - V2D-T- ~ W2D-T- «2D = V 2D -^— 

Re oy oz) oy 



(21) 



for the 0(1) streaks field t/2D- For a well-chosen z phase, this 
2D-3C field has X and S2 symmetries. 

The three-dimensional part of the initial guess was then de- 
termined by computing a neutral instability mode of the rolls 
and streaks field for the chosen value of a using the linear sta- 
bility solver. For a = 0.577, this mode, which has ^2 sym- 
metry, is marginally stable when the forcing amplitude is A ^ 
4.4. Using the total rolls and streaks field and a small three- 
dimensional neutral mode perturbation as an initial guess, we 
finally succeeded in computing an exact three-dimensional non- 
linear steady forced solution using the Newton solver. In order to 
obtain a solution for zero forcing, we then only had to reduce the 
amplitude of the forcing term step-by-step down to zero ampli- 
tude and to follow the corresponding forced solution using the 
nonlinear continuation code. The bifurcation diagram with re- 
spect to the forcing amplitude A is shown in Fig.Q] In this plot, 
the y-averaged amplitude of the k x — a mode of the shearwise 
vorticity 77 

Amp 3D = % J" (j]e\p(-iax))dy (22) 

is used to measure the amount of three-dimensionality of the so- 
lution. As shown in this diagram, the whole procedure works 
because the bifurcation to three-dimensional solutions is sub- 
critical with respect to the forcing amplitude. Note also that A 
is effectively 0(1) at the bifurcation point. 

Before cyclonic continuation was attempted, the non- 
rotating plane Couette flow solution obtained previously was fi- 
nally continued with respect to the Reynolds number in order 
to obtain the full bifurcation diagram for the chosen aspect ra- 
tio. The Reynolds number of the saddle-node bifurcation point 
is Re = 127.7 (see Fig.©. 
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Fig. 1. Bifurcation diagram of forced nonlinear solutions in non- 
rotating plane Couette flow with (a,/3) = (0.577, 1.15) with re- 
spect to the forcing amplitude A, for Re = 150. Amp 3D , de- 
fined by Eq. ( l22l . is used to measure the amount of three- 
dimensionality in the full nonlinear solution. The resolution 
(N x , N y , N z ) = (32, 24, 32) was used to perform this calculation. 
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Fig. 2. Bifurcation diagram for non-rotating plane Couette flow 
with (a,fi) = (0.577, 1.15) with respect to Re. The nondimen- 
sional wall shear rate \d (U) /dy\(y — +1) is used as a measure of 
nonlinearity (for the laminar flow, \d (U) /dy\(+l) = 1). 



5.2. Continuation in the cyclonic regime 

The bifurcation diagram obtained in the previous paragraph was 
then used to compute solutions in the Rq > region by per- 
forming continuation with respect to Rn- Solutions for several 
Reynolds numbers were continued from the lower branch of 
Fig. [2] A turning point is found systematically at a critical ro- 
tation number depending on the chosen Reynolds number. The 
solution then evolves back to zero rotation number and reaches 
the upper branch of the non-rotating plane Couette flow bifurca- 
tion diagram. For the forementioned maximum affordable reso- 
lution, solutions could be obtained for Reynolds numbers up to 
500. Note that we suffered from a lack of numerical resolution 
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for Re = 500. The continuation curve is incomplete in that case 
and is given for illustrative purposes only, because the corre- 
sponding results are probably not very accurate. Unfortunately, 
it is not possible to increase the resolution of this experiment fur- 
ther by a significant amount because of both memory and CPU 
requirements. For instance, the LAPACK linear solver has algo- 
rithmic complexity 0{N 3 ), where N is the system size, so that 
doubling the resolution in even only one direction would be pro- 
hibitive (computing one point of the continuation curves using 
the present resolution already takes several hours on a 3 GHz 
processor). Branches of solutions obtained at various Reynolds 
numbers are plotted in Fig. [3] and a snapshot of isovorticity sur- 
faces and streamlines for the solution at the turning point at 
Re = 303 and Rq = 0.017 is shown in Fig. [4] Anticipating 
a comparison with results for anticyclonic rotation, we also 
present in Fig. [5] a bifu rcation diagram for Re = 303 similar 
to the one obtained by iNagatal dl990l) . showing both the 2D- 
3C Taylor vortices branch in the linearly unstable weakly an- 
ticyclonic regime and its connection with the three-dimensional 
wavy Taylor vortices penetrating into the cyclonic region. 

3.5- 
3.0 - 




1.5 -,' 



i.o L , , , i , , , , i , , , , i , , , , i , , J 

0.000 0.005 0.010 0.015 0.020 

Rotation number 

Fig. 3. Branches of three-dimensional nonlinear solutions ob- 
tained in the cyclonic region for Re = 150 (full line), Re = 303 
(dotted line) and Re = 500 (dashed line). A turning point corre- 
sponding to a critical rotation number is found for any Re. The 
Re = 500 calculation suffered from a lack of numerical resolu- 
tion. 

It is interesting to compare the critical Rn - Re curve ob- 
tained from our computations with already published experi- 
mental and numerical work. Figure [6] shows that the critical Re 
obtained for the exact coherent structures is smaller than the crit- 
ical Re for the onset of turbulence in cyclonic experiments for 
Re = 150 and Re = 303 (note again that the Re = 500 point in 
that plot is purely illustrative due to the lack of numerical res- 
olution). The critical Re is seen to increase with increasing cy- 
clonic rotation, similarly to experiments. A possible interpreta- 
tion for these observations is that the values obtained in our study 
for (a,P) parameters that minimize the global critical Reynolds 
number for the non-rotating plane Couette flow provide lower 
bounds on the transition Re in the cy clonic regime as wel l. Since 
iTillmark & Alfredssonl d!996l) and E esur & Longarettil d2005l) 
used a different aspect ratio in their experiments, larger Reynolds 
numbers may be required in their set-up to obtain turbulence. 
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Fig. 4. Positive (light brown) and negative (dark blue) isovortic- 
ity surfaces and randomly selected streamlines for the nonlinear 
solution located at Re = 303 and R n = 0.017. 
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Fig. 5. Bifurcation diagram for Taylor vortices (full line) 
showing their connection with three-dimensional rotating self- 
sustaining solutions (dotted line) at Re = 303. The dash- 
triple dotted curve shows the centrifugal stability limit, (a, /3) = 
(0.577,1.15). 

An alternative explanation is that the exact coherent structures 
computed here do not allow for a precise determination of the 
turbulence threshold in a real experiment and only give a rough 
estimate for the critical Re, because the actual transition process 
involves not only these solutions but all their possible instabili- 
ties, as mentioned previously. 

6. Anticyclonic rotation 

Several arguments regarding the existence of turbulence in 
Keplerian flows have been based on a possible symmetry with 
respect to Rq = -0.5. In this section, we report some continua- 
tion experiments in the anticyclonic region, which notably show 
that the self-sustaining process present in non-rotating or weaky 
cyclonic shear flows (Rq - 0) cannot be transposed into the 
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Fig. 6. Critical Reynolds number for exact three-dimensional 
nonlinear solutions for (a, ft) = (0.557, 1.15) as a function of 
the rotation number Rq (+), compared to t he experimental data 
obtained by iTillmark & Alfredssonl (119961) (o) for the onset of 
turbulence and the numerical data by Komminaho et alj d 19961) 
(a). The Re = 500 point only gives a qualitative indication of 
the actual critical value because of the lack of resolution for the 
corresponding numerical e xperiment. The numerical results of 
iLesur & Longa retti ( 200 5b for shearing box sim ulations match 
those obtained by ITillmark & Alfredssonl (17996 ) and are not re- 
produced here. 



Rn =: -1 regimes and that more complex nonlinear structures 
have to be searched for to find subcritical transition (if any) in 
this region. 

6.1. Finding a starting point for continuation 

Nonlinear continuation methods require a starting point, which 
is usually a bifurcation point. This bifurcation can have a physi- 
cal origin, when a linear instability like the centrifugal instabil- 
ity is present. In that case, continuation is straightforward. When 
there are no linear instabilities, a possibility is to include an extra 
artificial term in the original equations to trigger a linear insta- 
bility that would otherwise be absent, and to study the subse- 
quent nonlinear dyna mics when this artificial term is gradually 
removed. This is how jNagatal d!990l) proceeded in order to find 
three-dimensional solutions of non-rotating plane Couette flow. 
By adding a small amount of anticyclonic rotation to the origi- 
nal Couette flow, he caused the background flow to become cen- 
trifugally unstable and followed the 2D-3C Taylor vortices (ax- 
isymmetric structures with small amplitude streamwise rolls and 
strong streamwise velocity) into the nonlinear domain. In this 
regime, Taylor vortices become unstable to three-dimensional 
wavy Taylor vortices at a critical Reynolds number depending 
on the amoun t of ro tation of the system. The interesting find- 
ing of Nagata ( 1990) was that he was able to follow this three- 
dimensional solution down to zero rotation, thanks to a proper 
nonlinear feedback of the three-dimensional structures on the 
streamwise rolls originally generated by the centrifugal insta- 
bility. In the previous section, we proceeded in a similar way to 
find non-rotating solutions, but instead of adding a small amount 
of anticyclonic rotation originally, we used the same method as 
IWaleffd (12003) to sustain the rolls field, by adding an artificial 



forcing term in the streamwise yorticity equation. A third possi- 
bility, found by lClever & Bussel (Il992l) . is to cause the Couette 
flow to become unstable to convection by imposing a supercrit- 
ical temperature gradient in the shearwise direction in order to 
force the rolls field, and then perform a nonlinear continuation of 
the secondary instabilities of the convection rolls when the tem- 
perature gradient is suppressed. In spite of their a priori rather 
different characteristics, all these approaches lead to the same re- 
sult. The reason for this is linear transient growth via the lift-up 
effect. All these methods sustain weak streamwise rolls which 
give rise to a finite amplitude streaks velocity field, as a con- 
sequence of the spanwise redistribution of streamwise velocity. 
The instability of streaks, whatever the origin of streamwise rolls 
is, results fro m the presen ce of inflection points in the spanwise 
streak profile dJone 31981 . 

In order to investigate the nonlinear dynamics in the rotat- 
ing Keplerian regime, a proper starting point is required too, 
and a major problem is to find it. Unfortunately, we did not find 
a nonlinear connection between the non-rotating plane Couette 
flow solutions and the strongly anticyclonic region. As shown 
in Fig. [5] the three-dimensional lower branch of solutions that 
crosses the R& = line joins the two-dimensional Taylor vor- 
tex branch at some very small negative Rq, and we had difficul- 
ties in following the three-dimensional upper branch for Rq < 
-0.5 for the parameters of Fig. |5] (the solutions on this branch 
turn into nonlinear travelling waves (non-zero phase speed) at 
Rn ^ -0.216, similarly to what has been observed in experi- 
ments dAndereck et al.ll986tlNagatall990ln . Thus, another start- 
ing point must be found in order to study the region of large neg- 
ative Rq. It seemed difficult to start directly from the Keplerian 
regime, which is very far from any linear instability point, so 
our idea was to start from more convenient rotating regimes and 
to try to continue possible solutions towards Rq = -4/3. We 
experimented with two ideas for R& — -1 or close to minus one. 

The first working hypothesis, inspired bv lNagatal (1 1 9901) . was 
to start from values of Rq larger than but close to -1 in order to 
trigger a centrifugal instability, and to study the nonlinear be- 
haviour of secondary instabilities of the Taylor vortices in that 
regime, when Rq is made closer to -1 (centrifugal instability ap- 
proach). The second technique, similar to that used in Sect. [5] 
to obtain solutions for Rn = 0, was motivated by a very strong 
analogy and near-symmetry regarding linear theory between the 
Rn = and Rn — -1 regime s . From the strictly linear point of 
view, lAntkowiak & Branched (120061) have shown that the lift-up 
effect, which occurs for Rq = 0, has an equivalent for R& = -1, 
christened the "anti lift-up" effect. This equivalence can be seen 
by considering the linear equations describing the evolution of 
the streamwise velocity field u and of the streamwise vorticity a> 
(the rolls field) in plane Couette flow with system rotation: 



du 1 
— = -(R n + l)v+ — Am, 
at Re 

da) du 1 

ot dz Re 



(23) 
(24) 



For Rq — 0, a streamwise velocity field is generated by redis- 
tribution of shear by the streamwise rolls velocity field, leading 
to strong transient amplification and the formation of streaks. 
There is no source term in Eq. d24b and any streamwise vorticity 
perturbation must decay viscously in this linear framework. For 
Rn = -1, strong streamwise vorticity can be generated tran- 
siently by the Coriolis force acting on the streamwise veloc- 
ity. The streamwise velocity equation d23l has no source term 
in that case because the Coriolis force acting on the rolls ex- 
actly compensates the redistribution of the mean shear by the 



F. Rincon, G. I. Ogilvie and C. Cossu: Self-sustaining processes in Rayleigh-stable rotating plane Couette flows 



9 



rolls, so that streamwise velocity perturbations decay viscously. 
Transient growth observed for Rn — -1 can also be interpreted 
as an epicyclic oscillation of infi nite period, since the epicycli c 
frequency is zero in this regime dAntkowiak & Brancherll2006l) . 
Pushing the analogy with the non-rotating case further and now 
considering nonlinear effects, the generation of finite-amplitude 
rolls by the anti lift-up effect may then allow for the development 
of three-dimensional instabilities that could behave in a simi- 
lar way to their Rn = counterparts and feedback on a weak 
streamwise velocity field, thus leading to a self-sustaining pro- 
cess for Rq = - 1 . The procedure in that case is to set directly 
Rn = -1, to artificially force a weak streaks field in order to am- 
plify streamwise rolls, and to follow the nonlinear development 
of the instabilities of the rolls when removing progressively the 
forcing term (forcing approach). 

6.2. Centrifugal instability approach 

We start by describing results obtained in the first working hy- 
pothesis. Using both linear stability codes and nonlinear con- 
tinuation, we identified streamwise-independent Taylor vortices 
bifurcating from the spanwise-rotating background Couette flow 
at various Rn larger than but very close to -1. In practice, the 
bifurcation to Taylor vortices occurs when the Taylor number 



Ta = -16Re 2 # n (l + Rn) 



(25) 



exceeds 1707, so that the critical Rn depends only slightly on 
Re at large Re. We chose /3 = 1.5585, which is the critical 
wavenumber for this instability (in the Taylor-Couette experi- 
ment, this would be j6 = 3.117 because 2d is used as unit of 
length in that case). The structure of the full velocity field of 
the Taylor vortices in this rotating regime is quite different from 
the one at small Rn- close to Rn = -1, the streaks field has a 
lower amplitude than for small Rn, which is a consequence of 
the absence of the lift-up effect. The secondary instabilities of 
the Taylor vortices close to the Rayleigh line Rn = - 1 are there- 
fore different from those at small Rn, whose inflectional nature 
is due to the shape and amplitude of the streamwise velocity field 
in that regime. This marks a fundamental difference between the 
two regimes and already hints that the argument on symmetries 
of solutions with respect to Rn = -0.5 is limited. Due to the Z, 
symmetry of the Taylor vortices, two steady three-dimensional 
instabilities can be found, with either £2 (wavy Taylor vor- 
tices, hereinafter WTV) or X. (twisted Taylor vortices, here- 
in after TTV) sym metry. These instabilities have been identified 
by^ 
by 



Nagatald 19861) and TT Vs have been observed experimentally 
Andereck et al. (11983b for Rn ~ -0.8. Finally, some steady 



subharmonic instabilities with spanwise wavelength twice the 
wavelength of the Taylor vortices and a n instability w ith non- 
zero phase speed are also present (see iNagatal 1 1 988b . These 
latter instabilities are poor candidates for continuation into the 
Rayleigh-stable regime because they hav e been reported experi - 
mentally for -0.5 > Rn > -0.8 only (see lAndereck et aljfl983l) . 
Consequently, they have not been considered in this study. 

Close to the Rayleigh line, the first unstable three- 
dimensional mode (when Re is increased) is the wavy Taylor 
vortex mode for the numerically accessible range of Re (see 
Fig. [7] for a comparison between the growth rate of WTVs and 
TTVs for some of the parameters used here). In Fig. [8] we 
show the corresponding bifurcation diagrams in the vicinity of 
R n = -1 for Re = 500 and Re = 1500. The bifurcations to 
2D-3C Taylor vortices are supercritical with respect to Rn and 
the associated branches of nonlinear solutions do not cross the 



Rayleigh line Rn = - 1 ■ Independently of Re, the secondary bi- 
furcations to three-dimensional WTVs are unfortunately also su- 
percritical and the corresponding branches of solutions do not 
approach the centrifugal stability limits either, unlike their coun- 
terpart at Rn ~ 0, in spite of having the same £2 symmetry 
(compare with Fig. [5). The same remarks apply to TTVs. This 
demonstrates clearly that the phenomenology of self-sustaining 
processes at Rq = cannot be transposed directly into the 
Rn = — 1 regime, even though analogous linear transient am- 
plification mechanisms exist in each of these regimes. 
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Fig. 7. Growth rate of the three-dimensional instability modes of 
the Taylor vortices as a function of the streamwise wavenumber 
a for R n = -0.97875, Re = 500 and f3 = 1.5585. Full line: 
wavy vortex mode. Dotted line: twisted vortex mode. Nonlinear 
continuation from the bifurcation point is initialized by adding a 
small amount of neutral modes with a = 1 to the Taylor vortex 
flow. 



6.3. Pitfalls of the forcing approach 

As mentioned earlier, using a forcing approach similar to that 
presented in Sect. appeared as a plausible alternative to 
the centrifugal instability approach in order to identify a self- 
sustaining process at Rn = -1. Such an attempt wa s fur- 
ther encouraged by the work of IWedin & KersweiH d2004l) . who 
managed to identify a self-sustaining process in non-rotating 
pipe Poi seuille flow by con t inuing forced nonlinear solutions, 
whereas iBarnes & Kersw ell ( 2000) had previously been unable 
to continue nonlinear solutions obtained for a rotating pipe 
Poiseuille flow down to zero rotation. 

The simplest (and possibly naive) way of tackling the prob- 
lem was to include an artificial forcing term 



A ( 2 it 
F x (y, z) = ^ W + — I cos/fe cos 



(f) 



(26) 



in the streamwise velocity equation in the Rayleigh-neutral Rn = 
-1 regime in order to sustain the following streamwise velocity 
field against viscous dissipation: 



A lny\ 

= ^ cos ' 6zcos (yj- 



(27) 
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Fig. 8. Bifurcation diagrams in the Rp_ 1 region for (a,/3) = 

(1, 1.5585) and Reynolds numbers Re = 500 (lower full, dotted 
and dashed lines) and Re = 1500 (upper full, dotted and dashed 
lines). The resolution is (N y , N z ) = (24, 48) for the Taylor vortex 
branches (full lines) and (N x , N y , N z ) = (16, 24, 48) for the WTV 
branches (dotted lines) and TTV branches (dashed lines). The 
dashed-triple dotted line is the Rayleigh line. Continuation was 
stopped at smaller values of for Re = 1500 due to insufficient 
numerical resolution. 

The forcing amplitud e A is assumed to be 0(1), as in the non- 
rotating problem (see Waleffe 19951 120031) . in order to generate 
possibly unstable finite-amplitude 0(1) streamwise thanks to the 
Coriolis force acting on the 0(1 /Re) streamwise velocity field. 
Note that unlike for Rq = 0, there is no need to solve numeri- 
cally a Stokes problem for the streamwise velocity here, because 
the boundary conditions on u make it possible to use exact ex- 
pressions like Eq. d26| ) for the forcing term. 

Initializing the Newton solver with the streamwise velocity 
field d27| i and streamwise forcing term ( f26b plus a streamwise 
rolls field with streamwise vorticity a> being the solution of the 
linear problem 

Aw = -Re— , (28) 
oz 

nonlinear steady 2D-3C forced solutions could be obtained for 
various forcing amplitudes (see Fig. |9). However, their proper- 
ties turned out to be quite different from the prediction of linear 
theory that the streamwise rolls velocity field be (9(Re) larger 
than the streamwise velocity field. One of the reasons for this 
difference is that the rolls feed back on the streamwise velocity 
field which generates them via the anti lift-up effect and sweep 
it close to the boundaries. In comparison, the corresponding ef- 
fect at Rq = would be the advection of axisymmetric struc- 
tures like rolls by the streaks field (the ud x operator), which 
is strictly zero due to the streamwise-independence of the rolls 
and streaks pattern. Therefore, the analogy between Rq — and 
Rn = — 1 breaks down when nonlinearity is included, at least 
for the forcing term used here (the precise reasons for this break 
down and the choice of a suitable forcing term will be further 
discussed in Sect. |7|. The problem becomes even more appar- 
ent when computing the stability of these 2D-3C solutions to 
three-dimensional infinitesimal perturbations with wavenumber 
a. For Re = 500 and (a,/3) = (0.8, 1.5585) for instance, bifurca- 
tion to three-dimensional solutions with Z, symmetry occurs for 



A = <9(10 3 ) only, whereas bifurcation occurs for A = 0(1) in the 
Rq -0 case. Furthermore, the first three-dimensional instability 
found for Rq = -1 is supercritical with respect to the forcing 
amplitude (Fig. [9]) and subsequent three-dimensional nonlinear 
solutions cannot be tracked down to zero forcing (compare with 
Fig.CDforfln = 0). 

This second numerical attempt confirms the conclusions 
drawn from the centrifugal instability approach, namely that a 
self-sustaining process close to the Rayleigh line, if any, must 
act quite differently from its non-rotating counterpart. It could 
be argued that our values of Re are too small to observe any- 
thing interesting in the anticyclonic regime, or that the chosen 
aspect ratio does not allow for three-dimensional self-sustained 
solutions. We performed the same kind of computations with dif- 
ferent aspect ratios for Reynolds numbers up to approximately 
1500 but did not notice any qualitative change in the bifurcation 
diagrams, which always show bifurcation to three-dimensional 
solutions at large A and supercritical behaviour with respect to 
the forcing amplitude. 




Forcing amplitude A 

Fig. 9. Bifurcation diagram with respect to the forcing amplitude 
A for/fo = -1, (a,P) = (0.8, 1.5585) andRe = 500. The full line 
curve is the branch of 2D-3C forced solutions and the dotted 
line curve is the three-dimensional branch of solutions with X. 
symmetry. The bifurcation to these three-dimensional solutions 
occurs for A = 834.7 and is supercritical with respect to the 
forcing amplitude. 



7. Asymptotic description of self-sustaining 
processes 

The numerical results presented in the previous section illus- 
trate the difficulty of identifying self-sustaining processes in 
Rayleigh-stable regimes. However, there is currently no mathe- 
matical proof that such solutions do not exist, so that it worth 
attempting to understand why they cannot be found in plane 
Couette flow on the Rayleigh line and under which circum- 
stances they could be found. In this section, we provide some 
analytical arguments based on an asymptotic description of the 
problem at large Re which, together with the numerical results, 
highlight the fundamental differences between the R& = and 
Rn = - 1 regimes in plane Couette flow. Starting with the non- 
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rotating case, we obtain a two-dimensional formulation for the 
solutions of the lower branch of the non-rotating self-sustaining 
process at very large Re. We then consider several possible 
asymptotic descriptions of three-dimensional solutions close to 
the Rayleigh line and demonstrate why the actual Taylor vor- 
tices generated in this regime in plane Couette flow cannot lead 
to a rotating self-sustaining process, unlike their —Rq «: 1 coun- 
terparts. We finally show that self-sustaining solutions on the 
Rayleigh-line may only exist if the 2D-3C part of the flow satis- 
fies nontrivial solvability conditions. 

7.1. Non-rotating Couette flow 
7.1.1. Basic equations 

Following Sect. [3] perturbations u = (u, v, w) and p to the non- 
rotating plane Couette flow Ub — y e x satisfy the nonlinear equa- 
tion 

(d t + yd x + u • V) u + ve x = -Vp + eAu, (29) 
together with suitable boundary conditions and the constraint 



V ■ u = 0. 



(30) 



The parameter e = 1 /Re will be used in the following to perform 
an asymptotic expansion. 

7.1.2. Streamwise averaging 

Let an overbar denote averaging in the streamwise (x) direction, 
for which periodic boundary conditions apply. Then u and p can 
be separated into "mean" and "wave" parts, i. e. 

u(x,y,z,i) = u(y,z,t) + u'(x,y,z,t), (31) 

p(x, y, z, t) = p(y, z, t) + p'(x, y, z, t), (32) 

with it' = and p' - 0. We also separate u into toroidal (stream- 
wise) and poloidal (shearwise and spanwise) parts, parallel and 
perpendicular to the averaging direction: 

u = ue x + «p. (33) 

The mean (toroidal and poloidal) and wave parts of Eq. d29l ) are 
then 

(d, + Hp ■ V) S + v = eAu + F x , (34) 
(d, + Hp ■ V) ii p = -Vp + eAu v + F p , (35) 
(d t +yd x +u-V)u'+u'-Vu+u'-Vu'+FW e x = -Vp'+eAu', (36) 
with 

V-h p = V-h'=0, (37) 

where 

F 



-u' ■ Vh' (38) 
is the mean force associated with the wave. We also define 

G = V x F. (39) 

Since H p is a solenoidal two-dimensional velocity field we can 
write it in terms of a streamfunction if/(y, z, t): 



Hp = Vx(>//e x ). 
The associated streamwise vorticity is 
co = -Ail/. 



(40) 



(41) 



We therefore replace equations ( 1341 ) and d35l l, taking the curl of 
the latter to eliminate p, with 



d(ip, u) 

OfU — — — - + v = eAu + F x , 



d t co ■ 



d(y,z) 
d(i/f, co) 



eAu + G x . 



(42) 



(43) 



d(y,z) 

The reason that streamwise-independent solutions cannot be 
self-sustaining is that the streamwise vorticity co satisfies an 
advection-diffusion equation with no source term (since G origi- 
nates from the stream wise-dependent wave parts of the solution). 
For a self-sustaining solution, the mean and wave parts must in- 
teract cooperatively. 

7.1.3. Asymptotic reduction 

Remembering that the streamwise velocity field is strongly am- 
plified by the lift-up effect and is consequently far larger than 
the poloidal velocity field, we seek steady asymptotic solutions 
in the limit of very large Reynolds numbers (e <k 1), of the form 



u = u (y,z)e x + dii(y,z) + ■ 

h' = eu[(x,y,z) H , 

p = ep\(x,y,z) + ■■■ . 
At first order, these are required to satisfy 

So) a . ._ 
37 — r- + <Wi = Au 

(y + u )d x u\ + (u\ ■ VS + v\)e x = 
At second order, we have 



-Vp' 



d(y,z) 



Acoi + G X 2, 



with 



o>\ — —Atffi, 
Vu[= 0, 



-V x h' ■ Vh' . 



(44) 
(45) 
(46) 



(47) 
(48) 

(49) 

(50) 
(51) 
(52) 



These equations are closed. Equation (|48| | for the wave is lin- 
ear and can be solved by Fourier analysis in x. All the equations 
are therefore reduced to a two-dimensional form. This approach 
leads to a simplified set of equations which are just a reduced 
version of the full problem of determining exact self-sustaining 
solutions. Since the wave involves a single Fourier harmonic, 
in some sense we have achieved a radical spectral truncation 
of the three-dimensional problem, yet this is fully supported by 
the asymptotic analysis. The two main advantages of this set 
of equations are that (i) the equations are two-dimensional; (ii) 
there are no small parameters and therefore no small scales to 
resolve, since the Reynolds number has been scaled out. The 
amplitude of the wave is arbitrary since it satisfies a homoge- 
neous linear equation. For a domain of fixed size, the amplitude 
of the wave needs to be chosen so that So has the right ampli- 
tude to give rise to a marginally stable wave. Alternatively, the 
amplitude of the wave can be left as a free parameter and the 
streamwise wavenumber (or periodicity length) adjusted to al- 
low the marginally stable wave to exist. 

The agreement between this analysis and the self-sustaining 
process in non-rotating plane Couette flow is demonstrated in 
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Fig. [TOland Fig.Qj] Figure [TO] shows the evolution with respect 
to Re of the ratio between the vertical (z) kinetic energy con- 
tained in the (k x , k z ) = (a, /?), (2a, p), (3a, /3) Fourier modes and 
in the (k x ,k z ) = (a,/3) mode at y = 0.707, for the non-rotating 
self-sustaining solutions of the lower branch of Fig. [2] At larger 
Re, the energy contained in the k x > a modes becomes neg- 
ligible compared to the energy contained in the fundamental 
k x — a mode, which confirms the prediction made in the pre- 
vious paragraph that the problem can be solved using only one 
Fourier harmonic asymptotically. Figure QT| shows that the ex- 
pansion d44b is very well satisfied by the non-rotating solutions 
of the lower branch. The streamwise velocity amplitude remains 
approximately constant and 0(1) for increasing Re, while the 
poloidal velocity amplitude decreases like 1 /Re, as expected for 
a first order quantity in e. 



a 



10" 



lO" 4 



10" 



10" 



100 



1000 



Kcvnolds number 



Fig. 10. Vertical (z) kinetic energy ratio between either (k x , k z ) = 
(a,j3) (full line), (k x ,k z ) = (2a, j3) (dotted line) or (k x ,k z ) = 
(3a, ft) (dashed line) modes and the (k x ,k z ) = (a,/3) mode 
(y = 0.707) for the lower branch of non-rotating self-sustaining 
solutions, as a function of Re. The bump obtained for Re > 500 
and k x = 3a can be explained by the fact that an error tolerance 
of 10" 6 for the energy was allowed in these computations and 
that the (N v , N z ) = (24, 32) resolution starts to be insufficient for 
such values of Re. 

One strategy for solving the reduced system would be sim- 
ilar to Waleffe's approach to the full problem. One can initially 
set the wave amplitude to zero and replace the forcing term G& 
by an artificial force of suitable form, then solve the coupled 
equations for the "roll" ifii and "streak" uq resulting from this 
forcing and adjust the wavenumber a to locate a marginally sta- 
ble wave. If the "feedback" G X 2 has the same sign and a roughly 
similar spatial form to the artificial forcing, then the solution can 
be continued by adding a small amount of wave and reducing the 
artificial forcing down to zero. If a is fixed in this process, the 
amplitude of the wave will be adjusted in a self-consistent way. 
However, we have not carried out such a procedure. 

7.2. Asymptotics on the Rayleigh line 

As demonstrated in Fig. [10] and Fig. [TT] the asymptotics for 
Rn — are very successful in describing how the actual the self- 
sustaining process operates in non-rotating plane Couette flow. 



a 



10" 



10"* I 

100 



1000 



Reynolds number 



Fig. 11. Maximum of streamwise velocity perturbations u (full 
line) and of poloidal velocity perturbations Vv 2 + w 2 (dotted 
line) for the lower branch of non-rotating self-sustaining solu- 
tions, as a function of Re. The dashed line corresponds to a func- 
tion proportional to 1 /Re. 



Recalling that the non-rotating solutions are connected to sec- 
ondary instabilities of Taylor vortices for small unstable rotation 
(Fig- ID, it is clear that the asymptotic expansion (l44b is satisfied 
for —Rn « 1 as well with the Coriolis force playing the role of 
F, provided that -Rn ~ e 2 (or equivalently that the Taylor num- 
ber defined in Eq. d25l l is fixed). This shows that considering the 
slightly unstable rotating regime allows to capture the essential 
properties of the non-rotating self-sustaining process. 

Here, we similarly analyse the properties of the actual Taylor 
vortex flow close to the Rayleigh line and its relations with solu- 
tions at Rn = -1, and demonstrate why continuation to Rp_ < -1 
regimes through three-dimensional solutions is not possible in 
the same way as in the Rn ~ region in plane Couette flow. 
Introducing a critical amount of rotation Rn = — 1 to study self- 
sustaining behaviour on the Rayleigh line, our equations become 



(d t + yd x + u • V) h + ue y = -Vp + eAu, 



(53) 



V ■ u = (54) 
(note the u e y term here instead of the ve x term in Eq. d29l)). Then 
(d, + u p -V)u- eAu + F x , (55) 
(d, + u p -V)u p + u e y = -Vp + eAUp + F p , (56) 



(d,+yd x +u-V) u' +u' Vu+u' Vm' -u' ■ Vu'+u' e y --Vp'+eAu' 

(57) 

with 

(58) 



V ■ Hp = V ■ u = 0. 



Writing u p = V x (ij/e x ) and oi = -Ai/r, this time we have 

o(y,z) 

d/ w — jr~ ? ~ d z u = eAoj + G x . (60) 
o(y,z) 

The reason that streamwise-independent solutions cannot be 
self-sustaining is now that the streamwise velocity satisfies an 
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advection-diffusion equation with no source term. Recalling that 
we are looking for a self-sustaining process in the Rn = -1 
regime, this time we would like to take advantage of the linear 
anti lift-up effect present in this regime to perform the asymp- 
totic expansion. The anti lift-up effect predicts generation of 
streamwise rolls 0(1 /e) larger than the stream wise velocity field, 
which led us to consider two possibilities for asymptotic expan- 
sions. 

7.2.1. 0(e) streamwise rolls 

A natural asymptotic expansion on the Rayleigh line is obtained 
by considering the symmetry of 2D-3C Taylor vortices with re- 
spect to Rq = -0.5. If (u,v, w) describes a Taylor vortex so- 
lution at Rq, it can be shown that (-Rn/(Rn + l)u,v, w) is a 
Taylor vortex solution at R' n = -(Rq + 1). For a given Taylor 
number, -Rn(Rn + 1) ~ e 2 , which means that if u is the stream- 
wise velocity component of a Taylor vortex at -Rq «: 1 , then 
the streamwise velocity component of the corresponding Taylor 
vortex at R' a = -(Rn + 1) close to -1 is proportional to e 2 u. 
For slightly supercritical Taylor numbers and -Rn <K 1, v and 
w are 0(e) and u is 0(1). In the symmetry transformation, the 
rolls field remains 0(e), while the u field becomes 0(e 2 ) close to 
Rn — - 1 . Put differently, for slightly supercritical Taylor num- 
bers, the streamwise velocity for —Rn ^ 1 is 0(1 /e) larger than 
the poloidal velocity, which is consistent with the lift-up effect 
present in the limit Rp_ = 0, while for Rn slightly larger than -1, 
the poloidal velocity field is 0(1/ e) larger than the streamwise 
velocity, which is consistent with the anti lift-up effect present 
in the limit Rn = - 1 . Applying this symmetry transformation to 
the Rq — asymptotic expansion d44l i. we obtain the following 
expansion at Rn = - 1 : 



ewpifj.z) + e 2 u 2 (y,z) + 



Since the amplitude of this 2D-3C flow is small compared to the 
amplitude of the background flow, no secondary instabilities can 
be expected. In fact, any marginal three-dimensional secondary 
solution (u', //) must obey 



yd x u' + u'e y = —Vp' 



(62) 



at leading order in that case. The rolls field is not involved in this 
inviscid linear equation, whose only known solutions are simply 
2D-3C marginal Taylor vortices on the Rayleigh line that are not 
compatible with the requirement that u' = 0. 

The previous analysis therefore highlights one of the ma- 
jor differences between small Rn and Rn -1 in rotating 
plane Couette flow: for —Rn <s 1, secondary instabilities of the 
Taylor vortex flow occur for Taylor numbers close to the critical 
Taylor number for the onset of the centrifugal instability because 
the Taylor vortex flow is already dominated by 0(1) streaks in 
the slightly supercritical regime. For instance, the bifurcation to 
three-dimensional solutions in Fig. [5] occurs for a Taylor num- 
ber which is less than two times supercritical. For Rn — -1 and 
similar weakly supercritical Taylor numbers, there is no such 
secondary instability because the Taylor vortex velocity field is 
only 0(e) in that case. The secondary instabilities computed in 
Sect. [6] only appear at highly supercritical Taylor numbers. For 
instance, the Taylor number for which a secondary instability 
occurs for Re = 500 in Fig.[8]is ~ 200 times supercritical. The 
nature of these instabilities remains unclear. They may be related 
to the presence of strong velocity gradients in thin shear layers 
and boundary layers that form at large Taylor numbers, when 
nonlinearities (notably advection of the streamwise velocity by 



the rolls) become important (see Fig. [TZl ). Similar conclusions 
can be drawn when using the forcing term defined in Eq. d26| ) 
instead of the centrifugal approach, because the structure of the 
corresponding 2D-3C forced solutions looks very much like the 
structure of the nonlinear Taylor vortices for Rn close to -1. 



1.0 



>, 0.0 




Fig. 12. Streamwise velocity u (color scale) and poloidal (v, w) 
velocity perturbations (arrows) of the nonlinear Taylor vortex 
flow for p = 1.5585, R a = -0.995 and Re = 5000, correspond- 
ing to Ta = 1 166 Ta c . 



(61) 7.2.2. 0(1) streamwise rolls ? 



As shown previously, an essential requirement for non-rotating 
self-sustaining process and three-dimensional dynamics for 
-Rn « 1 is that secondary instabilities of 0(1) structures are 
present. In the light of these remarks, we finally attempt to de- 
rive an asymptotic expansion for Rn = - 1 taking advantage of 
the anti lift-up effect, that would allow for bifurcations of an 
0(1) x-indep endent velocity field. We seek steady solutions of 
the form (see lOgilvie & Lubowl (120061) for a similar expansion 
in a different context): 



u = H p oO>,z) + eui(y,z) + e 2 u 2 (y,z) + ••• , 
u' = eu\(x,y,z) + ■■■■, 
p' = ep\(x,y,z) + ■■■ . 
From zeroth to second order, these are required to satisfy 

d(ip , ojq) 



d(y,z) 

_ d(\liQ,ui) _ 
d(y,z) 

d(\po,u\) d(i//i,tL>o) 



= 0, 



d z u\ = Ao>o, 



d(y, z) d(y, z) 
(yd x + Hpo ■ V) it, + u\ ■ Vh p o + u\ e y = -Vp[, 

d(ifro,Uz) d(4f\,u\) 



d(y, z) d(y, z) 



Ami +F x2 , 



(63) 

(64) 
(65) 



(66) 

(67) 

(68) 
(69) 
(70) 
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with 

cu = -A^o, (71) 

a>i = -Atfru (72) 

V • u\ = 0, (73) 
F 2 = -u\ • Vh'j. (74) 

Note that a zeroth-order equation must be satisfied here be- 
cause of the presence of nonlinear self-interaction terms in the 
streamwise vorticity equation, contrary to the Rn = case 
where the zeroth-order terms were strictly vanishing thanks to 
the streamwise-indepence of the zeroth-order streaks. There is 
also no equivalent to Eq. d67| i in the non-rotating problem. The 
meaning of these two zeroth-order and first-order equations is 
that advection by the rolls, although it does not affect the overall 
energy budget, dominates over viscous and Coriolis terms. 

We now drop the subscript on i/tq and introduce a right- 
handed non-orthogonal coordinate system (x, r, iff) based on 
the leading-order streamlines, where t is the travel time mea- 
sured along a streamline from its intersection with some refer- 
ence curve. The Jacobian determinant for this transformation is 
8(y, z)/8(t, ifi) = 1. The preceding equations then become 

d T cj = 0, (75) 

d T u x = 0, (76) 

d T u x -{8^)8^- d z u x =<9 T (g^wo) + W^<V<4>), (77) 

8 T u 2 - (c^OcVSi = + cV(/%Si) + Fx2, (78) 

where g'i are the inverse metric coefficients. The unknown quan- 
tities co\ and «2 are eliminated by integrating with respect to t. 
Defining the circulation on the streamline labelled by iff to be 

rw = J g^dr, (79) 

where the integral is along the full length of the streamline, the 
following solvability conditions now need to be satisfied: 

U r ^) + I idzai)dT = °' (80) 
d / duA r 

^ra + r- dT =°- (81) 

They state that the diffusion of the x-independent stream- 
wise vorticity must be balanced by the Coriolis force acting 
on the streamwise mean velocity u\, while the diffusion of 
the jc-independent streamwise velocity must be balanced by a 
wave flux (or by a streamwise force in the absence of three- 
dimensional motions). Finally, the first order equations for the 
"wave" read 

(yd x + «po ■ V) u\ + u\ ■ Viipo + u\ e y = -Vp[ , (82) 

V ■ u\ = 0. (83) 

These equations are linear and can in principle be solved by 
Fourier analysis in x, as in the Rn = case. 

To summarize, the asymptotic scenario at Rn = -1 can be 
realized if i) a particular distribution of streamwise vorticity sat- 
isfying Eq. d66*l l can be generated thanks to an artificial forc- 
ing term satisfying the solvability conditions (fSOlt - (T8TT> ii) the 
three-dimensional instabilities of this forced flow have the abil- 
ity to replace the force thanks to their nonlinear feedback on the 
streamwise-independent part of the flow. In order to satisfy the 



first requirement, it may notably prove necessary to break the 
point reflection symmetry 

(y, z) — > (-y, L z /2 - z), (u, v, w, u) — » (-«, — v, —w, ai), (84) 

of Taylor vortices which causes the streamwise "forcing" term 
averaged over closed streamlines in Eq. ( IBTT i to vanish. In the 
light of the developments presented in this section, relaxing the 
constraint of the presence of rigid boundaries that are respon- 
sible for the formation of boundary layers may also help to 
make progress. The fact that transition for Rn < 1 has only 
been observed in shearing boxes so far seems to support this 
view. However, we have been unable to isolate such a flow in 
the course of our investigations. As mentioned in Sect. [3] stan- 
dard continuation methods are not appropriate to compute solu- 
tions in shearing boxes and only theory o r direct time-s t eppin g 
calculations such a s those pe rformed by lHawlev et al.l d 19991) . 
iLesur & Longarettil d2005l) or lShen et ail ( 120061) may be able to 
provide some further insight in that case. An alternative possibil- 
ity would be to consider a remapping of the shearwise coordinate 
to obtain an unbounded shear flow between -oo < y < +oo and 
to use nonlinear continuation techniques for such a flow. 

8. Discussion and conclusions 

In this paper, we have tried to build on the current understand- 
ing of transition in non-rotating shear flows in order to iso- 
late exact nonlinear solutions of the Navier-Stokes equations in 
Rayleigh-stable rotating shear flows and to identify a possible 
route to turbulence in such flows, including Keplerian flows. On 
the one hand, we have shown that nonlinear mechanisms act dif- 
ferently on the Rayleigh line, even though analogous transient 
linear growth processes exist in both types of flows. We have 
notably computed exact nonlinear solutions for no rotation and 
for cyclonic rotation but have not succeeded in isolating solu- 
tions for Rn — -1 and smaller using similar techniques. Using 
the asymptotic approach presented in Sect. [7] we have shed 
some light on the differences between the two regimes and have 
demonstrated quantitatively that one cannot transpose straight- 
forwardly the nonlinear stability properties of non-rotating shear 
flows and the phenomenology of non-rotating self-sustaining 
processes into the Rn — -1 regime. In particular, our results 
show that transition in the narrow gap Taylor-Couette experi- 
ment beyond the Rayleigh line, if any, cannot occur through a 
self-sustaining mechanism analogous to the one pertaining to 
non-rotating plane Couette flow. On the other hand, our findings 
do not rule out subcritical transition in Rayleigh-stable shear 
flows. Finding a way to satisfy the requirements of Eqs. d66l l- 
(f80b-(T8n> may help to discover self-sustaining solutions on the 
Rayleigh line, which would then offer an interesting starting 
point for a possible continuation towards Rn = -4/3. We have 
argued in Sect. [7] that the presence of walls may prevent such 
solutions in plane Couette flow with Rq = - 1 and that solutions 
with the right self-sustaining asymptotic properties may only be 
found in unbounded or periodic domains. 

An alternative to a streamwise rolls-based scenario is that 
the transit ion relies on the existence of shearwise rolls, as sug- 
gested by ILesur & Longarettil d2005l) . However, such structures 
are likely to be transient due to the ambient shear, and it seems 
difficult to identify a steady self-sustaining process involving 
them using continuation. This possibility leads us to finally come 
back to the Keplerian problem and to discuss the re cent re- 
sults of iBalbus & Hawlevl d2006l) and lShen etaTI d2006l) regard- 
ing the nonlinear outcome of the transient am plification of non- 
axisymmetric shearing waves in this regime. IBalbus & Hawlevl 
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d2006l) have notably questioned transition to turbulence via this 
process by showing that shearing waves (even three-dimensional 
ones) do not have any self-interactions and are therefore exact 
nonlinear solutions of the equations. We would like to emphasize 
that nonlinear interactions of transiently amplified structures are 
not a necessary condition for subcritical transition. Bypass tran - 
sition in shear flows was highlighted bv lTrefethen et al.l (119931) . 
who argued that this mechanism could trigger turbulence if ini- 
tial perturbations could reach finite amplitude and interact non- 
linearly. But a very simple model for the no nlinear interactions 
was used at that time, and Waleffd (1995b showe d (an argu- 
ment similar to the one by Balb us & Hawlevl d2006l) ). using the 
full Navier-Stokes equations, that the finite-amplitude streaks 
generated by the lift-up effect did not feedback on the stream- 
wise rolls directly, leaving the outcome of transient growth and 
the existence of self- sustaining solu tions very uncertain. The 
real breakthrough of IWaleffd d 19951) . however, was to notice 
that these spanwise modulated structures became unstable to 
three-dimensional disturbances when they reached finite am- 
plitudes, which allowed him to discover the full self-sustaining 
process of which t his ins tabi lity is a crucial pie ce. In fact, both 
IBalbus & Hawlevl d2006l) and lShen et all d2006l) have found that 
the transient amplification of shearing waves leads to fairly com- 
plex dynamical behaviour as well, even though the nonlinear 
self-interactions of these structures are vanishing, because such 
waves become unstable to Kelvin-Helmholtz instabilities when 
they reach finite amplitude. These authors did not observe subse- 
quent self-sustaining behaviour for the range of parameters ex- 
plored in their simulations, but more investigations on this prob- 
lem may be needed before definite negative conclusions can be 
made. To be complete, it should also be mentioned that the pres- 
ence of a nonlinear saturation mechanism for transiently ampli- 
fied structures is not a sufficient condition for subcritical tran- 
sition either: one of the reasons why a simple self-sustaining 
process based on the anti lift-up effect cannot occur for wall- 
bounded shear flows in the Rq = -1 regime is precisely the 
nonlinear advection of the streamwise velocity field by linearly 
amplified streamwise rolls. Finally, we emphasize that follow- 
ing the nonlinear evolution of shearing waves does not ensure 
that the whole nonlinear dynamical behaviour in Rayleigh-stable 
regimes has been explored. For instance, the nonlinear cyclonic 
solutions reported here do not come out of the transient ampli- 
fication of non-axisymmetric shearing waves. Instead, they are 
connected to non-rotating solutions which owe their existence to 
the transient growth of axisymmetric structures. Even though we 
have not been able to isolate nonlinear coherent structures in the 
anticyclonic Rayleigh-neutral regime in this study, solutions in 
the Keplerian regime, if they exist, may well have a connection 
with nonlinear solutions at Rq = -1, for which linear axisym- 
metric transient growth should be p resent through the anti lift- 
up effect. The numerical results of iLesur & LongarettH d2005l) 
seem to indicate that transition at Rq = -1 is possible (at least 
in shearing boxes) and that a connection with Rayleigh-stable 
regimes exists, at least for a limited interval of Rq. 

Another debated point which has connections with our re- 
sults is that transition in disks could not be explained like transi- 
tion in wall-bounded shear flows because boundary layers play 
an important role in transition in the latter case, whereas they 
do not exist in accretion disks (or in shearing box simulations). 
IBalbus & Hawlevl d2006l) consider the case of subcritical travel- 
ling Tollmien-Schlichting waves in non-rotating plane Poiseuille 
flow with no-slip boundaries to illustrate this idea. Such an ar- 
gument would be important if transition in wall-bounded shear 
flows was known to be triggered by instabilities relying on the 



presence of walls. However, as noted in Sect. |2j transition in 
many shear flows (including plane Poiseuille flow with stress- 
free boundaries) cannot be explained this way simply because 
there is no such instability in these flows. The self-sustaining 
process is the most natural explanation available so far for tran- 
sition in non-rotating shear flows, but boundary layer s are clearly 
not an essential part of it. In fact, [Waleffel d2003l) has found 
that self-sustaining solutions exist for both no-slip and s tress- 
free boundary conditions and ILesur & LongarettiTd2005h have 
shown that transition in non-rotating or cyclonic shearing boxes 
and plane Couette flow with walls occur under very similar con- 
ditions. Also, note that the instability of the streaks field in the 
non-rotating case does not result from the shearwise inflection of 
the mean velocity profile (which is due to the presence of bound- 
aries at y = ±1), but has its roots in the spanwise (vertical in 
the disk terminology) inflections of the streamwise-independent 
part of the flow instead. This kind of inflection is not precluded 
in flows that are not wall-bounded. Of course, we have demon- 
strated in this study that what happens for Rayleigh-stable an- 
ticyclonic regimes is a completely different story, but the fact 
that transition for Z?q < -1 has only been observed for shear- 
ing box configurations and the results of Sect.|7]indicate that the 
presence of walls at least does not favour subcritical transition in 
these regimes. 

A lot of issues regarding subcritical transition to turbulence 
in rotating shear flows remain very unclear. There is currently 
no mathematical proof that there is no such route to tu rbu- 
lence for Rq < — 1 , and th e results of lHawlev et alj d 19991) and 
ILesur & LongarettH d2005l) indicate that turbulence indeed exists 
in this regime in shearing box configurations. If this is true, some 
structures must be present in the phase space of the correspond- 
ing nonlinear dynamical system in order to ignite the transition 
process, and it should be possible to identify them. We hope that 
the results reported in this paper, which have been obtained by 
combining linear and nonlinear approaches, will be helpful to 
achieve this goal. 
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